Chapter 19 

Operators and Observables of the Dirac field 



from my book: 
Understanding Relativistic Quantum Field Theory 

Hans de Vries 
November 10, 2008 



Chapter 



Contents 



19 Operators and Observables of the Dirac field 1 

19.1 The components of the Dirac field 2 

19.2 The total mass/charge operator 3 

19.3 The general form of the operators 3 

19.4 Space-like integration and Special Relativity 4 

19.5 The total mass/charge current operator 5 

19.6 The Haniiltonian time derivative operator 5 

19.7 The average position operator 7 

19.8 The average velocity operator 7 

19.9 The acceleration operator 8 

19.10 The intrinsic spin operator 13 

19.11 The orbital angular momentum operator 14 

19.12 Conservation of total angular momentum 14 

19.13 Plane wave angular momentum conservation 15 



Chapter 19 

Operators and Observables of 
the Dirac field 



2 



Chapter 19. Operators and Observables of the Dirac field 



19.1 The components of the Dirac field 

Having studied the propagation mechanism of the Dirac field in the pre- 
vious chapters using the chiral representation we can separate the most 
general form of the Dirac field in the following main components. 

J 

Where Q{t, x) is the scalar real valued amplitude field and |Qp determines 
the density, the exponent exp( i(f){t^ x) ) defines the local phase of the wave, 
and the bi-spinor {^l, ^r) describes the chiral components. The bi-spinor 
is normalized locally as. 



iR )\iR ) Uw Ui? 



CRiL + CiiR = 2mc2 (19.2) 



Note that the overline swaps the left and richt chiral components besides 
taking the complex conjugate values. From the same local bi-spinor we 
can obtain the local energy/momentum or charge/current density 



iR P \ Tr ^ ^ CR<y''^R + CL~^^iL = 2( E, p^c, pyc, p,c ) (19.3) 



with af" = ( a\ (T^, ), and a'' = ( a\ -a^, -a^, -a^ ) 



The above values are known locally at each point so that we can derive 
a local physical velocity {vx,Vy,Vz) for each point since Vi = c^pi/E. In- 
ternal consistency requires that the local phase can be expanded for an 
infinitesimal environment like. 

^icf>{t,x) _ ^-iEt/h + ip-x/h + ia (19 4) 



Which provides a second way to determine the four-momentum locally at 
each point by means of the phase change rates, (which can be determined 



19.2 The total mass/charge operator 
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independently from the magnitude change rates). Both the spinor compo- 
nents as well as the phase change rates must transform consistently going 
from one reference frame to another. The four-momentum derived from 
both of them must be equal. 



19.2 The total mass/charge operator 

This is the simplest operator possible. The mass and charge are Lorentz 
scalars. To obtain a Lorentz scalar via an integration over space we need 
compensate for the Lorentz contraction. The bilinear covariant which does 
this is the mass/charge-density ^p^^ip. 



m 



^ J dx' V^tV (19.5) 



qtot = ^ I dx' V'tV (19.6) 



The operators are merely constants which adjust the normalization. The 
mass/charge-density 'ip^^ip is generally normalized as 2mc^, using the rest- 
mass energy. 



19.3 The general form of the operators 

More involved operators act on the wave-function ip. They can act on the 
bi-spinor ^ using 7-matrices or on the phase/amplitude component using 
differential operators. If O is the operator then the general form is. 

O = / dx^ ^p-f^O^ (19.7) 



From the notation used above we can recognize that the operator acts on 
a Dirac field. We can also use another equivalent notation which resembles 
the way operators are applied in the non-relativistic Scrodinger theory. 



^ Cnorm 



J dx^ tp*dij (19.8) 
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The equivalence becomes clear when we write this out in more detail. The 
7^^ undoes the swap between the left and right chiral components of "ip 



19.4 Space-like integration and Special Relativity 

In the examples following in this chapter we will use the general form given 
below, where the O is a "local density" operator on the field. 

O = [ dx^ (19.10) 



Local because ip*Oip yields a result for each point and density because it 
acts on the "Lorentz compressed" field: is larger (denser) by a factor 
7 as i/jip which compensates for the Lorentz contraction of the total field 
during integration. 

There is a fundamental problem however with the spatial integration with 
which the reader will be familiar with from the version of expression in 
non-relativistic Quantum Mechanics. 



The issue is that physical space-like integration is forbidden in Special 
Relativity: It propagates information instantaneously from x = — oo to 
X = +00. We need to be aware of the limitations of the above expression. 



If the field represents a static time-independent wave-function then there 
is no physical objection to the integration. We can in general expect the 
result of the integration to be the same in different reference frames (hyper- 
planes) . 

However, if there is interaction and the field changes then we can not 
expect the result of the integration to be the same in different reference 
frames. 



19.5 The total mass/charge current operator 
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19.5 The total mass/charge current operator 

The second simplest operator possible gives us the total electric current 
or the total mass-fiux. These represent the three spatial components of 
four-vectors which have the Lorentz invariant mass and charge as zeroth 
component. 

We simply use the current-density ?/i7*'0 which forms a four- vector together 
with 'tpj^tp. 

^mass = ^Jdx^ V^TV (19.11) 

/ = [ dx^ ipYij (19.12) 

zmc^ J 

We find the operators O if we rewrite these in the general form 'ijj^^Otp 
discussed above. 

ijjf^ = ijjry^j'^j'tP = lP*dtp (19.13) 

& = = ( 7' ° ) (19.14) 



19.6 The Hamiltonian time derivative operator 

We define the Hamiltonian operator as the first order time-derivative in 
order to be able to use it as a commutator to obtain the time derivatives 
of other observables. We recall the Chiral 4d Dirac equation. 



/ 
/ 



dt 



+ c 





-a' 



a" 




d_ 



mc 
ih 



i^R 



(19.15) 



So we can define the Hamiltonian operator by multiplying all terms with 
the factor icHj^. 
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d . ^ f -a' \ d , / I ^ 2 



We can write the total Hamiltonian in an integral form as follows. 

H = jd.^ n = ^,jd.^ r ^ (19.17) 

Where H is the Hamiltonian density. Note that the time-derivative only 
acts on the exponential component of ^. It is sandwiched between '^^^'^ 
where it acts as the Hamiltonian-density operator. 

It acts as a density operator because the bilinear covariant ijj'^f^ijj transforms 
as the charge density. Note that the Lorentz transform of the charge 
density is expressed by the spinor product ^7°^, where iIj=^ e-xp{iijj) . It 
is the relativistic spinor transform rather than the transform of the phase 
change rates which determines the transform of the density. The density 
has to increase by a factor of 7 due to the Lorentz contraction. 

If is a solution of the Dirac equation then we can express the Hamiltonian 
density as both sides of (19.16) 



(19.18) 

We should be able now to obtain the time derivative of an operator O by 
means of a commutation with the Hamiltonian operator H as follows. 



dt h L 



H , O 



' (^HO - OH^ (19.19) 



Where the Hamiltonian operator H is given by. 

H {=ihdt) = -ich'j^-r'di + 7Vc^ (19.20) 



19.7 The average position operator 
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19.7 The average position operator 

The average position operator can be defined as well by using the mass/charge 
density field tpj^tp. 



The mass / charge density increases as 7 with the velocity because the area 
which confines the wave-function decreases with 7 due to the Lorentz con- 
traction. Both effects cancel when the density is integrated over space. 
The total mass/charge is Lorentz invariant. 



Here we have sandwiched x as an operator between ■07'^'0 as a "position- 
density" operator, which simply means that it yields the average position 
when integrated over space. We can place x anywhere since the three com- 
ponents of X act as constants in the (three) expressions for the components 
of ?/'7°x ip. 



19.8 The average velocity operator 

Having obtained a position operator, we are now able to derive a velocity 
operator by commutating it with the Hamiltonian operator. Recalling the 
position and Hamiltonian operators. 

X' = x' (19.23) 

The commutator of these two becomes a constant matrix for each of the 
three components. 



V V 



H,X' 
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The mass term did commute but the differential term didn't. Consistent 
with the other operators we now write for the average velocity. 

This is just, (see equation (19.14)) 

Where f is the current density. The current density increases by an extra 
factor 7 depending on the speed because the volume containing the moving 
electron becomes smaller by a factor 7 due to the Lorentz contraction. The 
integral over space of the current density gives the total current as we did 
see with the current operator in section (19.5). The total current J* is 
related to the total charge Q by the relation 

/ = VavgQ (19.28) 



demonstrating that (19.25) represents the correct velocity operator. 



19.9 The acceleration operator 

Having shown that the velocity operator yields the correct results when 
applied properly we can derive an acceleration operator with the use of the 
Hamiltonian. 



A' = 



~dt 



H, V 



(19.29) 



The velocity operator anti commutes which each term of the Hamiltonian 
with the exception of the momentum in the same direction. The commu- 
tation with the mass term of the Hamiltonian gives. 



19.9 The acceleration operator 
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/ 
/ 



mc , c 



— (T 



imc / cr' 
^ V -a^ 



(19.30) 



While the commutation with the momentum derivative terms of the Hamil- 
tonian results in. 



-ich 



-a^ \ d 








a'' \ 


d 


) 





a' 



(19.31) 



In more verbose notation we can write for . 



A' 



I 

h L 



imc 



(19.32) 



Where S'^' is the 4x4 matrix with two cr^ on its main diagonal. Sandwiching 
this between the field gives us. 



imc" 



(19.33) 



We should expect that ip*A^jp is zero for a plane wave so we will study 
the above expression with this in mind. We find that the above expression 
transforms like the Lorentz Force. 



m dt 



iq 
m 



(vxB + (19.34) 



With the E component zero in the rest frame and B pointing in the same 
direction as the spin in the rest-frame. We therefor associate E and B with 
the magnetization and polarization of the field itself due to the magnetic 
moment. 



m dt m \ €0 



(19.35) 
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The polarization of the field is the result of the Lorentz transform of the 

magnetic moment if the velocity is not equal to zero. This Lorentz force is 
zero in the case of the plane-wave. It the particle is at rest then only the 
magnetization is non-zero, there is no Lorentz force. 

The Lorentz Force evidently stays zero after transformation. The forces 
from the magnetization and the polarization cancel. In other words: 
A plane wave doesn't experiences acceleration. 

We now return our attention to the result obtained sofar. 

3 

tP*A'ip = ci^ (7S^ X 5i) ^ + V;7°7V (19.36) 

First we recognize ipj^'j^ip as the three spatial components of the axial 
vector current of the Dirac electron. 

j\ = V^tS'V' (19-37) 



Secondly we see that i/jdjip corresponds with the momentum of the plane 
wave. It acts only on the phase change rates of the exponential since the 
density and the bi-spinor are considered constants throughout the plane- 
wave. 

Thirdly we recognize 'tp'y^^^ip from the Gordon decomposition of the vector 
current as the polarization P* of the electron field. The polarization is 
the result of the magnetic moment of a boosted electron. The magnetic 
moment is Lorentz transformed into an electric dipole moment. 

- « = - ^ K "o' 

The Gordon decomposition involved the matrix a^'^ which contains the 
matrix-products used above in the zeroth row and column. 



19.9 The acceleration operator 
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a 








a 






-a' 











1(7 













-ay 





-1(7 









1(7 

















icr?' 
ia 



(19.40) 



It was shown that the bihncar tpa^^ip can be associated with the mag- 
netization/ polarization tensor of the electron's field due to the magnetic 
moment. 



2mc^ 






- Px 


- Py 


- Pz 









icMy 


Py 







-icM^ 


Pz 


-\CMy 








(19.41) 



So, we have associated the second term at the right hand side of equation 
(19.36) with the polarization of the field due to the transformation of its 
own magnetic moment. 

3 

i^-i^i'i) is identical to ^ P (19.42) 
h egm 

We now need to associate the first term at the right hand side of equation 
(19.36) with vxM. 



c i/j (7^7^ X di) ^ is identical to ^ [v x Mj (19.43) 
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In a plane- wave everything is constant except the phase cxp{—iEt/h+ipx). 
The differential operator act only on the phase so we can separate the two 
operators below. 

ii = V;7S'V P' = i'dj^ (19-44) 



The term is thus proportional to the cross product of and j\. What we 
need to show is that. 

v'xM^ cx p'x j\ (19.45) 



And this has to hold in all reference frames. The magnetization M* and 
the axial current density j\ point in the same direction as the spin in the 
rest-frame, however. The transform differently. The first transforms like as 
a combination of tensor component while j\ transforms as an axial vector. 

The transform of both NP and j\ depend on the relation of the pointer to 
the boost. The orthogonal and parallel components transform differently. 
Here we only need to take the orthogonal components into account due to 
the cross-product with the velocity/momentum, and indeed, the orthogo- 
nal component of M* transforms with a factor 7 more as the orthogonal 
component of j\ so, equation (19.45) holds. 



= = -(vxnoM+-p\ (19.46) 

m at m \ J 

In order to make the relation to hold we have used an extra condition here, 
namely that both the bi-spinor and the phases are identically Lorentz 
transformed. That they have both the same velocity. Without this as- 
sumption the plane wave would have a non-zero acceleration. 



19.10 The intrinsic spin operator 
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19.10 The intrinsic spin operator 

The operator for the intrinsic spin of the Dirac particle is given by. 

S' = hj^fj^ = (19.47) 



Where we can write out terms 7^7*7^ and expHcitly as. 



/ 

1 





-a' 







-/ 
/ 



M '"J ° ) (19.48) 



The matrix owns its name to the o"* matrices on the diagonal. The 
operator acts on the bi-spinor part of '0. Applying the operator gives 
us the total spin integrated over space. 



h 



s = 



2mc^ 



= — 



h 



2mc^ 



(19.49) 



The quantity j\='4>Y'y^ V' is the axial current density as we have seen in 
the chapter on the bilinear covariant fields of the Dirac particle. The spin 
and the axial current density j\ always point to the same direction, 
independent of the reference frame. 

The size of the components of increase with an extra factor of 7 because 
it is a density and the volume of the field decreases by a factor 7. If 
s represents the unit spin-vector at rest, and with s± and S|| being the 
components of s orthogonal and parallel to the boost then we can write 
the transform as. 



(s± +-5||7)a; 

{s± + S||7)2 



(19.50) 
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19.11 The orbital angular momentum operator 

The operator for the orbital angular momentum is analogous to the clas- 
sical version L = rx p and obtained by the substitution of with —ihdi. 



Q ( ydz- zdy 

zdx — xdz 
xdy - ydx 



V = -ihfx—^ = -ih\ zdx -xdz I (19.51) 



This operator acts on the phase/amplitude components of via the dif- 
ferential operator. The total orbital angular momentum integrated over 
space is given by. 

V = J dx^ ip*{-inrxdi)ij (19.52) 

The outer product with r can be brought outside the "sandwich". 

= 2S/ " ^ (^*^^^) ^^^-^^^ 



19.12 Conservation of total angular momentum 

When Dirac first published his theory in the beginning of 1928, [1],[2] he 
already had made an important discovery concerning the spin 1/2 nature 
of the electron. 



Studying the conservation of the orbital and spin angular momentum he 
found that the two were not independently conserved. He commutated 
the two operators with the Hamiltonian to obtain the time derivatives and 
found. 



dU _ i 
~dt ~ h 



dS' 
~dt 



h 





ich 7° 


f l^d. 


-Ydy 


H, V 


Ydx 


- I'^dz 








-l^dx 




= - 2ich 7° 


( l^d. 


-Ydy 


H, S' 


l^'dx 













(19.54) 



(19.55) 



19.13 Plane wave angular momentum conservation 
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Which we can write more condensed by using the cross product of the 7 
matrices and the partial derivatives. 

— = ichj^j'xdi), —= -2zc;i7°(y x5i) (19.56) 



Which shows that neither of the two is locahy conserved. The similar- 
ity between the two however leads directly to an expression for the total 
angular momentum operator which is locally conserved. 

r = D + (19.57) 



The time derivative of this operator obtained via the Hamilton is locally 
zero, so J* is a conserved quantity, or in other words " a constant of the 
motion" . 



(19.58) 



dt h 

This showed that the electron has an intrinsic magnetic moment oi \h 

19.13 Plane wave angular momentum conservation 

It turns out, on closer inspection, that in case of a plane-wave both spin 
and orbital angular momentum are in fact locally conserved independently. 



fifi f)C- 
— = ichj''{fxdi), —= ~2ichj'{j'xdi) (19.59) 



This is not directly obvious in the expressions above. We first have to let 
the operators act on the wave-function tp. 

ich ip* (7V xdi)ip = ich i) (7^ X di) il) (19.60) 
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The derivatives act purely on tlic phase-change rates of tlic exponential in 
the case of a plane-wave while the 7 matrices act only on the bi-spinor. 
We recognize two different currents in the expression above. First there is 
the Dirac vector field which represents the four-vector momentum of the 
particle. 

p' = ^V^tVtV = ^V^t'^V' (19-61) 

Secondly there is the expression which derives the momentum of the plane 
wave from the phase changes of the complex exponential. 

P' = (19-62) 

The momenta should be equal. Note that first is derived from the bi- 
spinor and the second from the phase-change-rates using the derivatives. 
Equation (19.60) represents the cross-product from the two which should 
be zero. 

Note that we used an extra condition here, namely that both the bi-spinor 
and the phases are identically Lorentz transformed. Although evidently a 
reasonable assumption, it is only with this extra condition that the spin 
and orbital angular momenta are conserved separately. 
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